CeCoIn 5 — a quantum critical superfluid? 



S. Ozcan, D. M. BrounJ^] B. Morgan, R. K. W. Haselwimmer, and J. R. Waldram 

Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge, CB3 OHE, United Kingdom 



O 
O 

(N 

a 

o 
o 

;h 

Or 

•4-3 



I 

-a 
a 
o 

o 



> 

On 

O 

VO 
O 
(N 
O 
-t— > 



i 

C 

O 

o 



J. L. Sarrao 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

Saeid Kamal, C. P. Bidinosti, and P. J. Turner 

Department of Physics and Astronomy, University of British Columbia, Vancouver B.C., 

(Dated: February 1, 2008) 



Canada V6T 1Z1 



We have made the first complete measurements of the London penetration depth A(T) of CeCoIns, 
a quantum-critical metal where superconductivity arises from a non-Fermi-liquid normal state. Us- 
ing a novel tunnel diode oscillator designed to avoid spurious contributions to A(T), we have es- 
tablished the existence of intrinsic and anomalous power-law behaviour at low temperature. A 
systematic analysis raises the possibility that the unusual observations are due to an extension of 
quantum criticality into the superconducting state. 



PACS numbers: 74.25.NF, 74.70.Tx 

A number of experiments on the recently discov- 
ered family of heavy fermion superconductors CeRhlns, 
Celrln5 and CeCoIns point to the existence of a non- 
Fermi-liquid (NFL) metallic state in these compounds. 
Central to this conclusion is the observation, in Celrlns 
and CeCoIns, that two of the key parameters of a Fermi 
liquid — the electronic heat capacity coefficient 7 = C/T 
and Pauli susceptibility \ — increase on cooling and show 
no sign of entering a temperature-independent Fermi- 
liquid regime |j| . Similar behaviour is observed in other 
heavy-fermion materials and is understood in terms 
of magnetic fluctuations near a zero-temperature criti- 
cal point, where theory predicts either 7 ~ -lnT or 
7 = 70 — AT 1 / 2 , depending on the dimensionality and 
nature of the magnetism || [I| . What happens to such a 
system when it becomes superconducting is at present an 
open question. Do the Fermi liquid parameters continue 
to evolve within the superconducting state, or does su- 
perconductivity abort the approach to quantum critical- 
ity? Also, would the absence of quantum-critical super- 
fluidity in materials such as the cuprate superconductors 
rule out a zero-temperature critical point as the source 
of their NFL normal-state behaviour? These are impor- 
tant issues, which we address in this Letter with the first 
complete and well controlled measurements of the Lon- 
don penetration depth A(T) of a superconductor known 
to be situated near a magnetic quantum critical point. 
Our novel oscillator design and carefully characterised 
samples avoid extrinsic contributions to the A(T) signal, 
and a systematic analysis leads us to consider new, NFL 
physics as the reason for the anomalous power-law be- 
haviour observed in X(T) in the low temperature limit. 

The normal phase of CeCoIns contains all the hall- 
marks of quantum criticality. De Haas-van Alphen 
(dHvA) measurements reveal large cyclotron masses that 
are strongly field dependent, exceeding 100 mo at low 



fields ||. The resistivity is reminiscent of the cuprates, 
nearly linear in temperature below 20 K ||. In mag- 
netic fields large enough to suppress superconductivity, 
C/T ~ — In T down to the lowest temperatures Jl|, |) . 
From 2.5 K up to 100 K the spin-lattice relaxation rate 
has the T 1 / 4 temperature dependence expected near an 
antiferromagnetic instability j?j , indicating proximity to 
a zero-temperature magnetic critical point. This is sup- 
ported by the fact that CeCoIns is a higher-density ana- 
logue of CeRhlns, a material in which weak, ambient- 
pressure antiferromagnetism gives way to non-s-wave su- 
perconductivity with T c = 2.1 K at a pressure of 1.6 GPa 

The superconducting state of CeCoIns is also highly 
unconventional and appears to bear strong similarities to 
that of the cuprate superconductors |l(| . Measurements 
of specific heat and thermal conductivity k(T) reveal low 
temperature power-law behaviour, similar to expecta- 
tions for a pairing state with line nodes || [ll| . A strong- 
power law temperature dependence is also observed in 
the low-temperature microwave surface impedance [ p"2| . 
The four-fold modulation of k by an angle-dependent 
basal-plane magnetic field indicates that the line nodes 
lie along the [110] directions ]T3]. This is confirmed by a 
field-rotation study of dHvA oscillations, which increase 
in amplitude below H C 2 for fields in the [110] direction 
In addition, observations of Pauli-limited supercon- 
ductivity in k(H ) p3] and a Knight shift that decreases 
below T c in all directions [m] imply spin-singlet super- 
conductivity. Together the experiments strongly suggest 
that pairing occurs in a d x 2_ y 2 state. 

CeCoIn 5 forms in a tetragonal crystal structure with 
alternating layers of Celn3 and CoIn2, and superconducts 
at ambient pressure below T c = 2.25 K. The high-quality 
CeCoIns crystals used in our experiment were grown by 
a self- flux method in excess In M H] . Although the start- 
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ing materials are very pure to begin with (Ce: 99.99%, 
Co: 99.9975% and In: 99.9995%), growth in excess In is 
expected to further refine the crystals. These naturally 
form as large a6-planc platelets, with mirror-like surfaces, 
and are ideally suited to high-frequency measurements. 
The high homogeneity and low defect level of the crystals 
are confirmed by our microwave measurements, which 
show a sharp superconducting transition (AT C < 30 inK) 
and a low quasiparticle scattering rate (1/r, the width of 
the conductivity spectrum ct(cj), rs 2 x 10 10 s _1 at 1.2 K) 
p2| . We focus on data from one sample, a crystal with 
dimensions a x b x c = 1.38 x 1.37 x 0.073 mm 3 . 

The A(T) measurements reported here were made with 
a 130 MHz tunnel diode oscillator (TDO) operated in 4 He 
and dilution refrigerator cryostats over the temperature 
range 0.1 K to 9 K. In addition, a thorough study of the 
microwave conductivity cr(u),T), obtained from surface 
impedance measurements from 1 to 75 GHz, was made 
on the same crystals used in the TDO experiments. Al- 
though the microwave measurements form an important 
part of this work (as a sample characterisation tool, for 
cross-checking and calibrating the TDO measurements, 
and as a means of determining the absolute penetration 
depth) a detailed discussion of those results will be pre- 
sented elsewhere 0|. We focus here instead on describ- 
ing the distinctive features of our tunnel diode oscillator 
and the results obtained using it. 

A good description of the general principles of oper- 
ation of the TDO is given in Ref. |lj, and oscillators 
of this sort have been used by other workers to study 
unconventional superconductors (l6[ 0. The chief in- 
novation in our apparatus is the novel geometry of the 
oscillator's probe circuit — a very small, self-resonant 
superconducting coil that acts as a local probe of the 
magnetic penetration depth. The miniature resonator 
has a square cross-section of side 0.5 mm and is wound 
from superconducting Nb wire on a high purity sapphire 
former. The high quality factor of the coil (Q — 2.5 x 10 5 
at 1.3 K) allows the resonator to be coupled inductively 
to the tunnel diode, with the mutual inductance backed- 
off until the diode only marginally sustains oscillation. 
This results in a frequency stability better than 1 part 
in 10 9 per hour and a field at the sample of ~ 10 nT. 
The combination of a high-Q resonator and low-power 
tunnel diode keeps the total heat load of the oscillator 
below 2 /j,W and allows the whole circuit to be cooled 
into the mK temperature range. In the experiments, a 
basal-plane face of the single-crystal sample is attached 
to one end of the sapphire former with vacuum grease, 
about 50 (im from the end of the coil. In this geometry, 
the coil locally induces currents that flow in a 0.5 mm 
wide ring in the centre of the crystal face, as shown in 
the inset of Fig. El This is important when studying the 
superfluid response of electrically anisotropic materials 
such as CeCoIns, for two reasons: first, only basal-plane 
currents flow, eliminating contamination of the signal by 
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FIG. 1: Temperature dependence of the penetration depth in 
CeCoIns. Inset: A calculation of the induced surface current 
density J(x,y) in the CeCoIns sample. Note that current 
flows in a ring, well within the boundaries of the crystal face. 

c-axis currents; secondly, this experimental geometry is 
inherently insensitive to nonlocal effects, as electrons in a 
quasi-2D metal are intrinsically confined to move parallel 
to the basal plane Q. 

In our setup there is thermal contact between the sam- 
ple and the resonator, with the result that the temper- 
ature of the entire oscillator is swept during the exper- 
iments. Being able to operate in this mode has the po- 
tential advantage that both resonator and sample can be 
embedded in a hydrostatic pressure medium, in princi- 
ple allowing the investigation of other quantum-critical 
systems, such as CePd2Si2 and Celn 3 [|l9), which only 
superconduct under high pressures. However, changing 
the temperature of the oscillator also introduces system- 
atic errors into the frequency shift signal. Fortunately, 
these are small over most of the temperature range and 
are highly reproducible, allowing them to be accurately 
accounted for using background measurements made in 
the absence of the sample. AA(T) is obtained from 
the oscillator frequency-shift signal A/o(T) = fo(T) — 
/o(Tbasc) using the cavity perturbation approximation, 
where AA(T) = -TAf (T). Here T = 6 A/Hz is 
a temperature-independent geometric factor determined 
empirically by comparison with penetration depth mea- 
surements made down to 1.2 K using a 1 GHz loop-gap 
resonator [ pp[ . 

The absolute value of A(T) has been determined from 
measurements of the surface impedance Z s = R s + iX s , 
made at 5.5 GHz with a TEon mode dielectric resonator. 
It is usually very difficult to determine A in this way, be- 
cause only shifts in X s with temperature are experimen- 
tally accessible. However, by carrying measurements of 
AX S (T) up to temperatures where the electronic scat- 
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FIG. 2: Two views of A(T) in the low-T limit, revealing a 
strong T dependence and the presence of gapless excitations. 
The solid line is a fit to the function AA(T) = AT 2 /(T + T*). 
Inset: A(T) is also well-described by a T power law. 



tering rate is much larger than the microwave frequency 
(i.e. for T > 20 K in CeCoIns) we access the Hagen- 
Rubens limit where R s ss X s . This enables the unknown 
offset in X s to be determined and, for sufficiently accu- 
rate AX S (T) data, gives X S (T) absolutely down to 1.2 K. 
(Confidence in our procedure is enhanced by the fact that 
R s and X s match between 20 K and 90 K, and that an 
analysis of the departure of R s from X s below 20 K reveals 
a temperature-dependence of the optical effective mass 
that follows the InT behaviour of C/T Q.) Having the 
X S (T) data is not enough: only in the low frequency 
limit, lot <C 1, does X s — io(1q\; at higher frequen- 
cies thermally excited quasiparticles also contribute to 
X s [|T). We use bolometric measurements of a(u>, T) ]2i| 
in the superconducting state to properly account for the 
quasiparticle contribution, obtaining A(1.2 K) = 3610 A 
and A = 2810 A in the T -> limit. This is a large 
penetration depth, characteristic of a metal with heavily 
renormalised electrons. 

Figure [l] shows the absolute A(T) data measured with 
the tunnel diode oscillator down to 0.1 K. In Fig. 0, a 
close-up view of the low-T data reveals a strong temper- 
ature dependence, indicating low-lying excitations and 
supporting the case for line nodes in the pairing state. 
However, A(T) does not have the simple linear T depen- 
dence expected for a d-wave superconductor and observed 
in the cuprates [^0) : instead the inset to Fig. || reveals 
that A(T) is better approximated by a T 15 power law, 
down to 0.1 K. (Similar curvature in A(T) has been in- 
ferred from 10 GHz X S (T) measurements and attributed 
to disorder Q . In those measurements, made down to 
only 0.25 K, lot ss 5 at low T and the clear connection 
between X s and the London penetration depth is lost.) 



FIG. 3: p B (T) = A 2 (0)/A 2 (T) for CeCoIn 5 , showing mean- 
field behaviour near T c and a strong, but nowhere linear, T 
dependence at low T. Inset: T dependence of the average 
effective mass m* required to bring the observed Ps{T) oc 
n s (T)/m* (T) into accord with the expected form for line 
nodes, n s (T) = no(l — aT/T c ), for several values of a. 



Also relevant is the observation that k(T) ~ 7 13 37 below 
0.2 K p| . The expected power law in the case of dilute 
strong-scattering impurities is n(T) ~ T 3 p3| , and an 
interesting possibility is that the fractional power laws 
in k(T) and A(T) might have the same origin. In Fig. ^ 
the data are plotted as p s (T) = A 2 (0)/A 2 (T). As p s is 
proportional to the ratio of the superfluid density n s to 
the effective mass m*, it is expected to have a linear T 
dependence in a d-wave superconductor, a geometric con- 
sequence of the presence of line nodes in the energy gap. 
This is not seen in our data, which show curvature over 
the full temperature range. 

We now consider possible explanations for the observed 
T-dependence of A, beginning with disorder. Strong- 
scattering impurities in a <i-wave superconductor are 
known to induce a crossover in A(T) from clean-limit, T- 
linear behaviour above T* w 0.56v / Wi7FT^ to a quadratic 
T dependence at low temperatures (Here m is the 

density of impurities and Tf the Fermi temperature.) We 
have assessed this possibility by fitting the interpolation 
formula of Ref. |§ AA(T) oc T 2 /(T + T*), to our data in 
Fig. |, and obtain T* = 0.3 K. Using T F = 50 K (from 
specific heat 0) we infer from T* that m = 0.26%. This 
is over an order of magnitude greater than the density 
of impurities in our starting materials. Thermal conduc- 
tivity experiments suggest the level of strong-scattering 
defects in our crystals is actually much lower: in k(T), 
T* < 30 mK, implying rij < 26 ppm [|ll| . 

Alternatively, Kosztin and Leggett have pointed 
out that nonlocal effects in a d-wave superconductor can 
cause an intrinsic crossover to T 2 behaviour in A(T), at 
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T* w 2T C £ /A , where £ is the BCS coherence length. 
Physically this is because the spatial extent of the Cooper 
pair is larger than Ao for a small range of angles around 
the nodal directions, and affects the field screening at 
correspondingly low temperatures. An estimate using 
the published value £ = 82 A §, gives T* rs 130 mK, 
significantly lower than the value obtained from the fit 
to our data in Fig. |[ In addition, for a quasi-2D d x 2_ y 2 
superconductor such as CeCoIn 5 , nonlocal effects should 
be very sensitive to geometry, being strongest for [100] 
surfaces, vanishing on [110] surfaces, and becoming ex- 
tremely small for [001] surfaces, where the 2D electronic 
structure forces the electrons to move almost parallel to 
the crystal face. Our penetration depth probe takes ad- 
vantage of this fact, only inducing currents on a [001] 
surface. In that case, T* is determined by the out- of - 
plane coherence length, £ c w 35 A ||, with the result 
that nonlocality should only become a consideration in 
our geometry below 56 mK. 

In a multi-band metal there is the possibility that su- 
perconductivity occurs strongly for only some pieces of 
the Fermi surface, and is induced in the others by an in- 
ternal proximity effect. This is thought to be relevant to 
CuO-chain superconductivity in the cuprates (2^, and 
to the physics of Sr2Ru04 p7| . In all cases it is expected 
to lead to a stronger T-dependence of A(T), due to the 
presence of a small temperature scale associated with in- 
trinsically weak superconductivity in parts of the Fermi 
surface. For CeCoIns, which appears to be a d x 2_ y 2 su- 
perconductor with line nodes, a proximity effect should 
result in positive curvature of p s (T) [p5| , opposite to that 
observed. 

The Bose-Einstein condensation theory of Ref. 



I pre- 
dicts AA(T) ~ T 1 - 5 at low T, and has been applied 
to A(T) measurements on the organic superconductor 
BEDT @. The theory has been developed for the 
case of nonoverlapping Cooper pairs in a short-coherence- 
lcngth, s-wave superconductor, as a possible explanation 
of the physics of the underdoped cuprates. It is not likely 
to be relevant to the case of CeCoIns, which has strongly 
overlapping, d-wave pairs. 

As none of the scenarios presented so far seems to pro- 
vide an adequate account of the A(T) data we propose an- 
other possibility, motivated by the NFL properties of the 
normal state. There we know that 7, x an d m * continue 
to be renormalised down to the lowest temperatures. If a 
T-dependent renormalisation also took place within the 
superconducting state, what would its effect on p s {T) be? 
This question is complicated by the fact that the differ- 
ent bands in CeCoIns have (field-dependent) cyclotron 
masses spanning the range 10 to 100 m S, and that 
p s oc (1/m*). Nevertheless, an effective mass that in- 
creased with decreasing temperature would introduce the 
right curvature into p s (T). To illustrate this point bet- 
ter we have taken our p s (T) data, assumed the d-wave 
form for the superfluid density, n s (T) = Uq(1 — aT/T c ), 



and plotted, below 0.8 K, in the inset of Fig. 3, the tem- 
perature dependence of m* for different values of a. In 
each case m*(T) shows a cusp-like upturn at low T, rem- 
iniscent of the quantum-critical behaviour of C /T in the 
normal state. Are there other data to support this con- 
jecture? A T-dependent renormalisation in the supercon- 
ducting state should show up most clearly in the specific 
heat. In CeCoIns the measurement is complicated by the 
large low-T entropy of the In nuclei but, when the nuclear 
contribution is subtracted from the data, C /T still devi- 
ates from the expected d-wave T-linear behaviour, even 
showing signs of a low temperature upturn p], [IlJ . The 
unusual low-T power law in the thermal conductivity, 
k(T) ~ T 3,37 , might also find a consistent explanation 
within this scenario. (As with p s (T), k(T) would be dom- 
inated by the least divergent band.) Further measure- 
ments of thermodynamic properties in the low-T limit 
should help assess the validity of this proposal. 

In conclusion, we have presented the first com- 
plete measurements of the London penetration depth of 
CeCoIns in the low-T limit and observed strong power- 
law behaviour that confirms the presence of low-lying 
excitations but is inconsistent with standard models of 
d-wave superconductivity. By working with carefully 
characterised samples and using a novel resonant-circuit 
probe that is inherently insensitive to nonlocal effects 
and electronic anisotropy we have been able to rule out 
extrinsic contributions to the A(T) signal. We propose 
an alternative interpretation, and suggest that the NFL 
renormalisation occurring in the normal state of CeCoIns 
might also take place within the superconducting phase, 
leaving us with the exciting possibility that CeCoIns may 
be the first example of a quantum-critical superfluid. 
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More recent measurements from E. Chia et al., which 
have come to our attention while preparing this work, 
confirm AA(T) ~ T 15 at low T, using crystals from the 
same source as ours. 
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